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Abstract 

A 1+1 dimensional hydro dynamical model in the light-cone coordinates is used to describe 



> 

o 

00 . central heavy-ion collisions at ultrarelativistic bombarding energies. Deviations from Bjorken's 



scaling are taken into account by choosing finite-size profiles for the initial energy density. The 
sensitivity of fluid dynamical evolution to the equation of state and the parameters of initial state is 
investigated. Experimental constraints on the total energy of produced particles are used to reduce 



\^ \ the number of model parameters. Spectra of secondary particles are calculated assuming that the 

O 

transition from the hydro dynamical stage to the collisionless expansion of matter occurs at a certain 



freeze-out temperature. An important role of resonances in the formation of observed hadronic 
O , spectra is demonstrated. The calculated rapidity distributions of pions, kaons and antiprotons 

in central Au+Au collisions at ^/s^ = 200 GeV are compared with experimental data of the 

•rH , 

^ ■ BRAHMS Collaboration. Parameters of the initial state are reconstructed for different choices 

of the equation of state. The best fit of these data is obtained for a soft equation of state and 
Gaussian-like initial profiles of the energy density, intermediate between the Landau and Bjorken 
limits. 

PACS numbers: 12.38.Mh, 24.10.Nz, 25.75.-q, 25.75.Nq 
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I. INTRODUCTION 



High-energy heavy-ion colhsions provide a unique tool for studying properties of hot 
and dense strongly-interacting matter in the laboratory. The theoretical description of such 
collisions is often done within the framework of a hydrodynamic approach. This approach 
opens the possibility to study the sensitivity of collision dynamics and secondary particle 
distributions to the equation of state (EOS) of the produced matter. The two most famous 
realizations of this approach, which differ by the initial conditions, have been proposed by 
La.da„ Q (fun .opp,,,) and B,o..e„ U (partial .ansparency). In .cent decade, 
versions of the hydrodynamic model were developed, ranging from simplified 



and 2+1 dimensional models 
to more sophisticated 3+1 dimensional 
mention the multi-fluid models {20], [21I 



mode 
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the Landau or Bjorken type 
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191 ]. One should also 



25l . l26l | which consider the whole collision 



process including the nuclear interpenetration stage. Recent theoretical investigations show 
that fluid-dynamical models give a very good description of many observables at the SPS 
and RHIC bombarding energies (see e.g. Ref. 2^]). 

The 2+1 dimensional hydrodynamical models have been successfully applied 
, [ll|, [12, [1^ to describe the px distributions of mesons and their elliptic flow at midra- 
pidity. These models assume a boost-invariant expansion [2] of matter in the longitudinal 
(beam) direction and, therefore, cannot explain experimental data in a broad rapidity region, 
where strong deviations from the scaling regime have been observed. More realistic 3+1 di- 
mensional fluid-dynamical simulations have been already performed for heavy-ion collisions 
at SPS and RHIC energies. But as a rule, the authors of these models do not study the 
sensitivity of results to the choice of initial and final (freeze-out) stages. On the other hand, 
it is not clear at present, which initial conditions. Landau-like [1^ or Bjorken-like j2|, are 
more appropriate for ultrarelativistic collisions. 

Our main goal in this paper is to see how well the fluid-dynamical approach describes the 
RHIC data on vr, K, p distributions over a broad rapidity interval, reported recently by the 



BRAHMS Collaboration 



28, 



291] . Within our approach we explicitly impose a constraint on 



the total energy of the produced particles which follows from these data. 

For our study we apply a simplified version of the hydrodynamical model, dealing only 
with the longitudinal dynamics of the fluid. This approach has as its limiting cases the 
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Landau and Bjorken models. We investigate the sensitivity of the hadron rapidity spectra 
to the fluid's equation of state, to the choice of initial state and freeze-out conditions. 
Modification of these spectra due to the feeding from resonance decays is also analyzed. 
Special attention is paid to possible manifestations of the deconfinement phase transition. 
In particular, we compare the dynamical evolution of the fluid with and without the phase 
transition. 

The paper is organized as follows: in Sect. [TTl we formulate the model and specify the 
equation of state, the initial and freeze-out conditions. In this section we also explain 
how we calculate particle spectra and take into account feeding from resonance decays. In 
Sect, mil we present the numerical results. We first consider the evolution of fluid-dynamical 
profiles for different equations of states and initial conditions. Then these results are used to 
calculate spectra of pions, kaons and antiprotons. In the end of this section we analyze the 
initial conditions motivated by the Landau and Bjorken models. Summary and outlook are 
given in Sect. HVl Some details of numerical scheme and the procedure to calculate feeding 
from resonance decays are given in Appendix. 

A short version of this paper is published in Ref . [s^ . 



II. FORMULATION OF THE MODEL 
A. Equations of ideal fluid dynamics 

Below we study the evolution of highly excited, and possibly deconfined, strongly- 
interacting matter produced in ultrarelativistic heavy-ion collisions. It is assumed that 
after a certain thermalization stage this evolution can be described by the relativistic fluid 



dynamics. Conservation laws o: 
following differential equations 



the baryon charge and 4-momentum are expressed by the 



|3l|: 



d,inUn = 0, (1) 

9,r'^" = , (2) 

where n, and T^^ are the rest-frame baryon number density, the collective 4-velocity 
and the energy-momentum tensor of the fluid. In the limit of small dissipation, one can 
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represent T^'^ in the standard form 54 1 



(3) 



Here e and P are the rest-frame energy density and pressure. 

We consider central colhsions of equal nuclei disregarding the effects of transverse col- 
lective expansion. In this case one can parametrize If^ in terms of the longitudinal flow 
rapidity Y as f/^ = (cosh F, 0, sinh F)^. It is convenient to make transition from the usual 
space-time coordinates t,z to the hyperbolic (light-cone) variables |2(], namely, the proper 
time r and the space-time rapidity rj , defined as 
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In these coordinates the equations ([I])-® take the following form |32 | 
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0, 
0, 



(4) 

(5) 
(6) 
(7) 



To solve Eqs. (l5])-(I7j), one needs to specify the EOS, P = P{n,e), and the initial profiles 
n(ro, 77), e(ro, 77), F(ro, 77) at a time r = tq when the fluid may be considered as thermody- 
namically equilibrated. 

In this paper we consider only the baryon-free matter, i.e. assume vanishing net baryon 
density n and chemical potential /i. In this case Eq. ([5]) is trivially satisfied and all thermo- 
dynamic quantities, e.g. pressure, temperature T and entropy density s = (e + P — nfi)/T , 
can be regarded as functions of e only. The numerical solution of Eqs. (I6l)-(1Z1) is obtained by 



using the relativistic version 



331] of the flux-corrected transport algorithm 3J]. The details 



of the computational procedure are given in Appendix A. 



B. Initial conditions 



Following Ref. 17|, we choose the initial conditions for a finite-size fluid, generalizing the 
Bjorken scaling conditions: 

(1^1 -^0)^ 



Y{To,r])=r], e(ro,?7) = eoexp 



2a2 



-Qi\v\-Vo) 



(8) 
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where 0(x) = (1 + sgnx)/2. All calculations are performed for tq = Ifm/c, assuming 
that Y and e are independent on the transverse coordinates. The particular choice 7^0 = 
corresponds to the pure Gaussian profile of the energy density. The case cr = corresponds 
to the "table-like" profile when the initial energy density equals eo at |?7| < r^o and vanishes 
at larger \ri\ . It is obvious that if a or tjq tends to infinity, one gets the limiting case of the 
Bjorken scaling solution. In this follows from Eqs. ©-©lX ~ ^ ^y^yP ^ire 

functions of r, determined by the ordinary differential equations 2] 55 1 



= , + P = . (9 

dr dr 

Comparison of initial energy density profiles corresponding to the same total energy of 
produced particles, but having different parameters a and r^o, is made below in Fig. [51 The 
symmetry of central collisions of equal nuclei with respect to the transformation z —z 
implies that at fixed r, including the initial state at r = tq , the flow rapidity and energy 
density should be, respectively, odd and even functions of t] : 

Y (r, -v) = -Y (r, r/) , e (r, -r]) = e (r, r]) . (10) 

Therefore, it is sufficient to solve Eqs. ©-(IT]) only in the forward "hemisphere" ?7 ^ , using 
the boundary conditions Y = 0, de/drj = at rj = . 

C. Equation of state 

One of the main goals of experiments on ultrarelativistic heavy-ion collisions is to study 
the deconfinement phase transition of strongly-interacting matter. In our calculations this 
phase transition is implemented through a bag-like EOS using the parametrization suggested 



in Ref. [l3[. This EOS consists of three parts, denoted below by indices H,M,Q, and 



corresponding, respectively, to the hadronic, "mixed" and quark-gluon phases. As already 



mentioned, pressure, temperature and sound velocity, Cg = a/ dP/de, of the baryon-free 
matter can be regarded as functions of e only. It is further assumed that Cg is constant in 
each phase and, therefore, P is a linear function of e with different slopes in the above- 
mentioned regions of energy density. 

The hadronic phase corresponds to the domain of low energy densities, e < en, and 
temperatures, T < ■ This phase consists of pions, kaons, baryon-antibaryon pairs and 
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hadronic resonances. Numerical calculations for the ideal gas of hadrons (see e.g. 35l |) 
predict a rather soft EOS: the corresponding sound velocity squared, = ~ 0.1 — 0.2 , 
is noticeably lower than 1/3. The mixed phase takes place at intermediate energy densities, 
from en to tq or at temperatures from Th to Tq . The quantity eq — eu can be interpreted as 
the "latent heat" of the deconfinement transition. To avoid numerical problems, we choose 
a small, but nonzero value of sound velocity cm in the mixed phase. The third, quark-gluon 
plasma region of the EOS corresponds to e > or T > Tg. It is assumed that = Cq 
reaches the asymptotic value (1/3) already at the beginning of the quark-gluon phase, i.e. 
at T ~ Tq . 

The analytic expressions for pressure and temperature as functions of e in all three phases 



are given by the following equations 



56] 
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P 



P 
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{eh <e< eq) 



(e > eg) • 



(11) 

(12) 

(13) 



Here Bm and Bq are the "bag" constants in the mixed and quark phases, respectively. Our 
parameters ch-, cm, cq are close to those used in Refs. ^, Q The constants Bm, Bq and Tq 
are found from the continuity conditions for P{e) and T(e) . The corresponding formulas for 



TABLE I: Parameters of EOSs with the deconfinement phase transition. 





en 

(GeV/fm3) 


eg 

(GeV/fm3) 








Th 
(MeV) 


(MeV) 


Bm 
(MeV/fm^) 


Bq 
(MeV/fm^) 


EOS- 


-I 


0.45 


1.65 


0.15 


0.02 


1/3 


165 


169 


-57.4 


344 


EOS- 


~II 


0.79 


2.90 


0.15 


0.02 


1/3 


190 


195 


-101 


605 


EOS 


III 


0.45 


1.65 


1/3 


0.02 


1/3 


165 


169 


-138 


282 



the entropy density are obtained from the thermodynamic relation s = {e + P)/T . Below 
we use several versions of this EOS with the parameters listed in Table [H 
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FIG. 1: Temperature dependence of energy density calculated from Eqs. (Ilip - (ll3p with pa- 
rameters from Table HI The solid line corresponds to both EOS-I (Tc = 167 MeV) and EOS-II 
(Tc = 192 MeV). The dashed-dotted line is calculated for EOS-III (T^ = 167 MeV). Open (full) 



dots show the lattice data 



3g] for baryon-free matter with number of quark flavors Nj = 2 (3) 



The parameters Th and Tq define the boundaries of a mixed phase region separating 
the hadronic and quark-gluon phases. The critical temperature as defined by lattice 
calculations should lie between Th and Tq , i.e. Tc ~ (Th+Tq)/2 . Earlier lattice calculations 
(see e.g. Ref. [37|) predicted the values = (170 ± 10) MeV for the baryon-free two-flavor 
QCD matter. However, a noticeably larger value = (192 ± 11) MeV has been reported 



recently in Ref. [38|. To probe sensitivity to the actual position of the phase transition, 
we consider the EOSs with different T^ and Tq (see Table [T]). The EOS-I corresponds to 
Th = 165 MeV and the parameters en, eg used in the parametrization LH12 of Ref. [l3| [57 1. 
In the EOS-II we choose Th = 190 MeV and scale en, eg to get the same values of e/T^ 
as a function of T /Th 58|]. The EOS-III differs from the EOS-I by choosing a higher value 



= 1/3. Finally, the parameters Bq,Tq are found from the continuity conditions for P 
and T . As one can see from Fig. [H the parametrization f|TT|) - f|T3|) of the EOS leads to the 
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FIG. 2: Comparison of different EOSs used in this paper. The hnes for EOS-I, EOS-II and 
EOS-III are calculated using Eqs. pT]) -(fT3|l with parameters from Table HI 

cusp in e/T^ at T = , typical for bag-like models. The lattice results for Nj = 3 are 
overestimated by about 30% in the temperature region from Tc to l.STc. An attempt to 
remove this discrepancy has been made in Ref. j^. It is also seen that the EOS-III strongly 
disagrees with the lattice data in the hadronic phase at T < Tc . 

Unless stated otherwise, these EOSs are used in the calculations presented in this paper. 
For comparison, we have performed also calculations for several purely hadronic EOSs. In 
this case we extend Eq. ffTTj) to energy densities e > with the same eu^Tu as in Tabled 
but taking different values of c\j from 0.15 to 1/3. In Figs. [2H3] we compare our EOSs 
with the phase transition and two hadronic EOSs with constant sound velocities Cg = ch ■ 
According to Figs. [H the mixed phase region in the EOS-II occupies a larger interval of 
energy densities, i.e. this EOS has a larger latent heat as compared to the EOS-I and 
EOS-III. By this reason, the life-time of the mixed phase should be longer for the EOS- 
II, assuming the same initial state. One can see that the hadronic EOS P = cj^e with 
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FIG. 3: Entropy density as a function of energy density for different EOSs. 

= 0.15 (the dashed hne in Fig. [2]) is much softer at high energy densities as compared 
to the EOS-I and EOS-IL However, according to Fig. [HI the difference between entropy 
densities predicted by this hadronic EOS and the EOSs with pase transition is rather small 
at e < 10 GeV/fm^. Due to these reasons, the rapidity spectra of pions and kaons predicted 
by our model for soft EOSs, corresponding to cj^ < 0.2, are, in general, rather insensitive 
to the deconfinement phase transition. On the other hand, the EOS of massless particles, 
P = e/3 , differs significantly from other EOSs shown in Fig. [3l As a consequence, visible 
differences with purely hadronic calculations appear only for hard EOSs, i.e. for cjj > 1/3 . 



D. Total energy and entropy of the fluid 

The equations of ideal fluid dynamics ([l])-(l3l) imply that the total baryon number B, 
total 4-momentum P^^ and entropy 5* are the same at any hypersurface lying above the 
initial hypersurface (in our case r = tq): 

B = j da^nU" = const , (14) 
9 



pM = y da^T^^ = const , (15) 
S = daiysU" = const . (16) 



In the case of central collision of identical nuclei, due to symmetry reasons, the total three- 
momentum is zero in the cm. frame. Therefore, Eqs. (flSj) give only one nontrivial constraint, 
namely, the conservation of the total cm. energy, E = . 

Let us consider a cylindrical volume of matter ("fireball") expanding in the longitudinal 
direction and choose the hypersurface of constant r . The elements of this hypersurface may 
be written as 

da^ = S^{dz, 0, dt)'' = T S^dr] {cosh T], 0, sinhr/)'', (17) 

where 5"^ = ttR^ is the transverse cross section of the fireball. Substituting (fT7|l into 
Eqs. f fT4|) - f|T6l) we get the relations 

+ 00 

B = S±T j dr]n cosh {Y -T]) , (18) 



-oo 

+ 00 



E = S_i_T J drj [e cosh Y cosh {Y — 7]) + P sinh Y sinh {Y — t])], (19) 



-oo 

-l-oo 



S = S^T jdris cosh (F - r/) . (20) 



for any fixed r ^ Tq . One can determine constants in Eqs. (fT4l) - ffT6l) by using Eqs. (fT8|) - (l20|) 
for r = tq. In the baryon-free case {n and B are zero), substituting r = Tq, Y = rj 
into ([I9])-(I20D leads to 

oo 

E = 2S±tq J dr] e{TQ,T]) cosh 1] , (21) 



oo 

S = 2S^ro [d7]s{To,v). (22) 







Equations ( j2T]) - (!22l) with E and S from (fT9l) - (|20l) can be considered as two sum rules for the 
fluid-dynamical quantities. We have checked that our numerical code conserves the total 
energy and entropy at any hypersurface r = const on the level better than 1% even for 
T > 200fm/c. 
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FIG. 4: Contours of constant eo in the (t — ijq plane, corresponding to the fixed total energy (I23p of 
initial fireball in central Au+Au collision {R = 6.5 fm). Points A-D represent the parameter sets 
included into the detailed analysis (see below). Point A (D) corresponds to the Gaussian (table- 
like) initial condition. Thin solid line shows the contour of total entropy S = 3 ■ 10^ , calculated 
for the EOS-I. 

Below we use Eq. (ETi) to constrain possible values of the parameters characterizing the 
initial state. This is possible since the total energy of produced particles can be estimated 
from experimental data. Indeed, the value of the total energy loss, AE = 73 ± 6 GeV per 
participating nucleon, has been obtained from the net baryon rapidity distribution in most 



281]. This gives an estimate of the total energy of secondaries in 



central Au+Au collisions 
the considered reaction: 

E = iVpart AE ~ 26.1 TeV , (23) 

where A/'part — 357 is the mean number of participating nucleons. 

Substituting the parametrization ([8]) into Eq. fl2T|) and taking the value of E from Eq. fl23l) . 
one gets the relation between the parameters eo, f]o, cr. The integration in the r.h.s. of Eq. (fT9|) 
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FIG. 5: Initial energy density profiles corresponding to the total energy E = 26.1 TeV in central 
Au+Au collisions. Solid and dashed-dotted lines correspond, respectively, to the parameter sets 
A and C (see Table IIIl|) . Dashed line is calculated for the table-like profile with r/o = 3.34, 
eo = 7GeV/fm^ (point D in Fig.|4|). Dotted line corresponds to the Gaussian profile with a = 1.42 , 
eo = 20 GeV/fm3 . 

can be performed analytically. This gives 

A = a((j) cosh ?7o + /?(ct) sinhr^o . (24) 

Here A = -E'/(2S'_LeoTo) and a,P are determined as 

a = y^o-exp ((tV2) , (3 = 1 + aeif (^-^^ . (25) 

It is easy to see that the solution of Eq. with respect to r/o exists at a < A , i.e. at not 
too large a. This solution is given by the expression 



Vo = In j . (26) 

Figure H] shows contours of constant eo in the tiq — a plane. All curves are calculated 

for the same total energy of the fluid given in Eq. (|23|) . The points at the horizontal 
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FIG. 6: Same as Fig. [5l but for initial temperature profiles. All curves correspond to the EOS-I. 

(vertical) axis correspond to Gaussian (table-like) initial energy density profiles. One can 
see that these profiles become narrower for higher maximal energy densities eo (in this case 
the contours of equal eo move closer to the origin of the a — tiq plane). It is obvious that 
an additional constraint on the total entropy S would further reduce the freedom in the 
choice of parameters eo, a . For example, the best fits of the BRAHMS data (see below) 
with the EOS-I fall on the line corresponding to 5 ~ 3 • 10^ . According to our analysis, 
the absolute yields of secondary mesons and antibaryons increase with even if the total 
energy is fixed. This effect is well-known from the Landau model [l| which postulates that 
the number of secondaries is proportional to S . The initial temperatures for the considered 
profiles reach the values 250-350 MeV, well above the critical temperature of the quark- 
gluon phase transition. 

Figure [5] represents some profiles of initial energy density used in this paper. All these 
profiles correspond to the same total energy of secondaries (l23l) . The corresponding profiles 
of the initial temperature, T(ro,?7), are shown in Fig. El They are calculated by using the 
EOS-I. Note, that at nonzero a the temperature profiles look much wider as compared to 
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those shown in Fig. [51 One can see that three different phases of matter appear aheady at 
the initial state. 



E. Particle spectra at freeze out 

The momentum spectra of secondary hadrons are calculated by applying the standard 



Cooper-Frye formula 39|, assuming that particles are emitted without further rescatter- 
ings from the elements da^ of the freeze-out hypersurface r = Tpiri) . Then, the invariant 
momentum distribution for each particle species is given by the expression 

where = (E,p)^ is the 4-momentum of the particle, y = tanh.~^(pz/ E) and px are, re- 
spectively, its longitudinal rapidity and transverse momentum, g denotes the particle's sta- 
tistical weight. The subscript F in the collective 4-velocity U , temperature T and chemical 



potential fi implies that these quantities are taken on the freeze-out hypersurface 591]. The 
plus or minus sign in the r.h.s. of Eq. (12 7p correspond to fermions or bosons, respectively. 

As has been already stated, the effects of transverse expansion are disregarded in our 
approach. Due to this reason, we cannot describe realistically the px spectra of produced 
hadrons, and analyze below only the rapidity spectra. For a cylindrical fireball with trans- 
verse cross section S± = ttR"^ expanding only in the longitudinal direction, one can write 
da^ = S± {dzp, 0, dtp)'^ . Using Eq. (jl]) we get the following relations 

d(Tf,P^ = S±{EdzF -Pzdtp) = S±mT {tf{v) cosh(?/ - 1]) - T'p{r]) sinh(?/ - r])} dr] . (28) 

Here rriT is the particle's transverse mass defined as rriT = a/ + p|, , where m is the 
corresponding vacuum mass. Substituting ( l28l) into (1271) and integrating over pt we get the 
expression for the particle rapidity distribution 



dN gSj 



dy 47r^ 



+ 00 

J dr] [rpir]) cosh(?/ — r]) — T'p{r]) sinh(?/ — i])] x 



X ldprpr,nr {exp (!!!1£2!!L^^Z_|M^) ± (29) 
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where Yp^r]) and Tpirj) are the flow rapidity and the temperature at r = tf^t]) . Note, that 
the Bjorken model 2] corresponds to Yp = rj and tf,Tf independent of rj. As can be seen 
from Eq. (12^ . the rapidity distributions of all particles should be flat in this case. 



F. Feeding from resonance decays 

In calculating particle spectra one should take into account not only directly produced 
particles but also feeding from resonance decays. In this paper we explicitly take into 
account the decays: p —>■ 27r, K* —>■ Kn, A Ntt , which are most important for describing 
the spectra of pions, kaons and antiprotons, respectively. The other resonances are taken 
into account using an approximate procedure explained below. The spectrum of i-th particle 



produced in the two particle decay channel R — > iX is calculated by using the expression 

oo 

^d^Nn^i d^Nn^i f , , , f . d^N^ Br^, JpPr ^ \ 

E—To = -7-0 — 1- = / dmRwimji) d PR— 61 Eoi ] , 

d-^p d-^prdy J J d'^pR Afiq^ \mR J 



40| 



(30) 



where the flrst integration corresponds to averaging over the mass spectrum of the resonance 
(mthr = 1^1 + mx), Pr = {Er,pr) and p = {E,p) are, respectively, the 4-momenta of the 
resonance R and the particle i . The coefficient bR^i denotes the branching ratio of the 
R iX decay. The energy and momentum of the i-th particle in the resonance rest frame 
are determined by the expressions 

,2 I ^2 ^2 



g.. = V.^? + ,g^ (31) 

The freeze-out momentum spectrum of the resonance, d^NR/ d^pR, is calculated using 
Eqs. (!27l) - (!28ll with m = rriR, g = qr. Below we assume that the freeze-out temperatures 
for directly produced particles and corresponding resonances are the same. Calculation 
of integrals in the r.h.s. of Eq. (1501) and the parametrization of w^uir) are discussed in 
Appendix B. 

The detailed calculations in this paper are made for spectra of charged pions {i = tt^), 
kaons {i = K^) and antiprotons {i = p). In these cases the most important contributions 
are given by decays of resonances R = Rq = p, K*{892), A (1232) , respectively. To take 
into account other hadronic resonances {R ^ Rq) we assume that the contribution of the 
resonance R is proportional to its equilibrium density, ur (Tp) multiplied by d}^, the average 
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TABLE II: Equilibrium densities of tt^-, iT^ -mesons and antiprotons hidden, respectively, in 
/9, K* and A resonances as well as corresponding enhancement factors at different temperatures Tp. 



Tf (MeV) 




(fm-3) 






* (fm-3) 


a^t 




(fm-3) 


P 


100 


5.74- 


10-4 


2.30 


1.08- 


10-4 


1.21 


1.48 • 


10-5 


1.09 


110 


1.37- 


10-3 


2.28 


2.85 • 


10-4 


1.26 


5.30 • 


10-5 


1.12 


120 


2.87- 


10-3 


2.31 


6.50- 


10-4 


1.30 


1.56 • 


10-4 


1.15 


130 


5.43 • 


10-3 


2.37 


1.32 • 


10-3 


1.35 


3.92 • 


10-4 


1.18 


140 


9.48 • 


10-3 


2.46 


2.46- 


10-3 


1.40 


8.74- 


10-4 


1.21 


150 


1.55 • 


10-2 


2.57 


4.24 • 


10-3 


1.46 


1.77- 


10-3 


1.24 


160 


2.41 • 


10-2 


2.70 


6.90- 


10-3 


1.51 


3.29 • 


10-3 


1.28 


165 


2.95 • 


10-2 


2.77 


8.63- 


10-3 


1.54 


4.39 • 


10-3 


1.29 
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number of i-th particles produced in R ^ iX decays (e.g. 
i = TT"*"). Within this approximation one has 



E 

R 



dyd'^PT 



a 



Ro—>i 



Ro 



dyd pt 



2/3, d' = 0.65 ... for 



(32) 



where the enhancement factor a^j^^ is 



a 



Ro 



\- d}^ nniTp) 



(33) 



Here ur is the total density of resonances R at the temperature T summed over all isospin 
states of R . For simplicity, we calculate in the zero width approximation. More details 
of riR and d^^ calculations can be found in Ref . 4l(| . We have checked for several resonances 
with two-body decays (e.g. for i = tt and R = K*) that such a procedure yields a very 
good accuracy. We include meson (baryon and antibaryon) resonances with masses up to 
1.3 (1.65) GeV and widths F < 150 MeV. The statistical weights, masses and branching 
ratios of these resonances are taken from Ref. 42] . The results for different Tp are shown 
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T (MeV) 

FIG. 7: Equilibrium densities of vr"*" mesons hidden in resonances rj,p,A,A as functions of tem- 
perature T . Solid line shows sum over all mesonic and baryonic resonances (see text). Dotted line 
represents the equilibrium density of free vr^-mesons (not hidden in resonances). 

in Table HTl One can see that factors a^j^ significantly exceed unity (especially for pions) and 
increase with temperature. 

Figure [7] shows the temperature dependence of vr "'"-meson densities, cf^ hr, hidden in 
different resonances R . One can see that, compared to other resonance decays, the channels 
p — s> 27r give the largest contribution at all temperatures. According to this calculation, 
in equilibrium hadronic matter the fraction of "free" pions (not hidden in resonances) is 
smaller than 50% at temperatures T > 150 MeV. 

III. RESULTS 

A. Dynamical evolution of matter 

Space-time evolution of matter as predicted by the present model is illustrated in 
Figs. [HHTSl In particular. Figs. [8H9] show profiles of the temperature and the collective 
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FIG. 8: Temperature profiles at different proper times r calculated for initial conditions corre- 
sponding to the parameter set A. Only the forward hemisphere (r/ ^ 0) is shown. Left and right 
panels correspond, respectively, to the EOS-I and the hadronic EOS P = cj^e with cjj = 0.15 . 

rapidity at different proper times r . Here we consider the Gaussian-like initial conditions, 
with parameters from the set A (see Fig. [Hand Table UTTl) . For comparison, the results are 
presented for the EOS-I and for the hadronic EOS with = 0.15. One can see that in the 
case of the phase transition the model predicts appearance of a a fiat shoulder in T{ri) and 
local minima in Y [r]) which are clearly visible at r < lOfm/c. This is a manifestation of 
the mixed phase which exists during the time interval Ar ~ lOfm/c. According to Fig. [HI 
the largest volume of this phase in the r^-space takes place at r ~ 5fm/c. In the considered 
case the "memory" of the quark phase is practically washed out at r > 30fm/c. As one 
can see from Fig. [9l at such late times deviation from the Bjorken scaling {Y = if) does not 
exceed 5%. 

Another important quantity which shows deviation from scaling hydrodynamics is the 
entropy density s multiplied by r. As noted in footnote to Eqs. ([9]), sr is constant in the 
Bjorken model. In Fig. [10] we show profiles of st at different r ^ tq , again for the parameter 
set A. Two conclusions can be made here. First, these profiles are only weakly sensitive to 
the presence of the phase transition. We have already noted this when discussing Fig. |3l 
Second, one can clearly see that st ^ const as a function of r . In contrast to the Bjorken 
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FIG. 9: Same as Fig. [51 but for collective rapidity profiles. 




FIG. 10: Same as Fig. [8l but for profiles of entropy density multiplied by r . 

model, ST drops by about 15% for \ri\ < 1 at t > 20fm/c. Due to this reason, we think 

which assume Bjorken 



that the above-mentioned 2+1 dimensional models 
scaling in the beam direction, are not very accurate even for the slice around rj = 0. 

To study sensitivity of flow to the EOS of the hadronic phase, we have performed calcu- 
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FIG. 11: Collective rapidity profiles at different proper times r calculated for initial conditions 
with parameters eo = 50 GeV/fm^, rjQ = 0, a = 0.98 . Left and right panel correspond, respectively, 
to the EOS-III and the hadronic EOS P = e/3. 

lations for different sound velocities ch keeping other parameters of the EOS the same as in 
Figs. [814T01 It is found that at fixed total energy E , acceleration of the fluid is stronger for 
harder EOSs characterized by larger ch ■ Figure [TT] shows profiles of flow rapidity for the 



EOSs with cjj = 1/3 60|. Comparison with Fig. [9] shows that the transition from = 0.15 
to = 1/3 results in about 3-fold increase of relative flow rapidity Y — rj (this quantity 
characterizes deviations from the inertial expansion of matter). Comparison of the left and 
right panels in Fig. [H] reveals that the sensitivity of flow to the phase transition is much 
stronger as compared with the soft EOS considered above {cjj = 0.15). The differences in 
acceleration at late times of 30-50 fm/c are especially visible for \f]\ ^2 . Qualitative expla- 
nation of this effect is obvious: fluid elements at such rj exhibit smaller acceleration in the 
case of phase transition. This is a manifestation of the mixed phase which is characterized 
by small gradients of pressure. 

Figure [I2] shows the fluid isotherms in the t] — r plane. The profiles of flow rapidity 
at hypersurfaces of constant temperature are shown in Fig. [131 This information is used 
to calculate particle momentum distributions by using formulae of Sect. IIIEI According 
to Fig. [12] (see the left panel), the initial stage of the evolution when matter is in the 
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FIG. 12: Isotherms in the r] — t plane calculated for the parameter set A. Left and right panels 
correspond to the same EOSs as in Fig. \8\ Shaded region indicates the mixed phase. 
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FIG. 13: Same as Fig. [121 but for profiles of collective rapidity corresponding to fixed temperatures. 

quark-gluon phase, lasts only for a very short time, of about 5 fm/c. The region of the 
mixed phase is crossed in less than 10 fm/c. This clearly shows that the slowing down of 
expansion associated with the "soft point" of the EOS plays no role, when the initial state 
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lies much higher in the energy density than the phase transition region. In this situation the 
system spends the longest time in the hadronic phase. The freeze-out at Tp = 130 MeV 
requires an expansion time of about 60 fm/c at 77 = 0. This is certainly a very long 
time which is seemingly in contradiction with some experimental findings. Indeed, the 
interferometric measurements 43| show much shorter times of hadron emission, of the order 
of 10 fm/c. As follows from our results, this discrepancy can not be removed by considering 
other EOS or initial conditions. A considerable reduction of the freeze-out times can be 
achieved by including the effects of transverse expansion and chemical nonequilibrium 
However, this will not change essentially the dynamics of the early stage (r < 10 fm/c) 
when expansion is predominantly one-dimensional. A more radical solution would be an 
explosive decomposition of the quark-gluon plasma, proposed in Ref . \^ . This may happen 
at very early times, right after crossing the critical temperature line, when the plasma 
pressure becomes very small or negative. We shall consider this possibility in a forthcoming 
publication. 



B. Rapidity spectra of secondary hadrons 



Below we show the results for rapidity distributions of vr- and fC-mesons as well as 



antiprotons produced in central Au+Au collisions at _jjs 



compared with data of the BRAHMS Collaboration 
lisions. 



28|, 



200 GeV. These results are 



29| for most central (0%-5%) col- 



TABLE III: Parameters of the initial states which give the best fits of the pion, kaon and antiproton 
rapidity spectra observed in central Au+Au collisions at ^Smn = 200 GeV. All sets correspond to 
the EOS-I. To denotes the maximum temperature at r = tq . Ei and £"3 are total energies of 
produced particles within the rapidity intervals |y| < 1 and \y\ < 3, respectively. 



set 


eo (GeV/fm3) 


a 


Vo 


To (MeV) 


El (TeV) 


E3 (TeV) 


E/S (GeV) 


A 


10 


1.74 





279 


1.53 


9.25 


0.89 


B 


9 


1.50 


0.62 


271 


1.54 


9.59 


0.86 


C 


8 


1.30 


1.14 


263 


1.49 


9.55 


0.86 
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FIG. 14: Rapidity distribution of vr"'" -mesons in central Au+Au collisions at ^/sJ^ = 200 GeV. 
Left panel shows results of hydrodynamical calculations for the EOS-I and initial conditions with 
parameters eo = lOGeV/fm'^, rjQ = 0, a = 1.74 (set A in Table Hn]) . Right panel corresponds to the 
EOS-II and the parameters eo = 5GeV/fm^, tjq = 0, a = 2.02. Solid, dashed and dashed-dotted 
curves correspond to different values of the freeze-out temperature Tp . The dotted line shows 
contribution of resonance decays in the case Tp = 130 MeV. Experimental data are taken from 
Ref. y. 



We have considered different profiles of the initial energy density, ranging from the 
Gaussian-like (t^q = 0) to the table-like {a = 0) initial states. We found that in the case 
of EOS-I it is not possible to reproduce the BRAHMS data on the pion and kaon rapidity 
spectra in Au+Au collisions by choosing either too small (eg < 5GeV/fm^) or too large 
(^0 ^ 15GeV/fm^) initial energy densities. For such eo values the pion and kaon yields can 
not be reproduced with any Tp . According to our calculations, the quality of fits is notice- 
ably reduced for initial energy density profiles with sharp edges, corresponding to a < 1 . 
As follows from the constraint fl25]) . such profiles should have a very large eo or a significant 
plateau — r/o < V < Vo- This would lead to more flat rapidity distributions of pions and 
kaons as compared to the BRAHMS data (see the next section). 

A few parameter sets which lead to good fits with the EOS-I are listed in Table IIIII 
In these calculations we choose various eo and a and determine rjo from the total energy 
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FIG. 15: Same as Fig. [T^ but for "'"-meson rapidity distributions. 

total entropy S* ~ 3 ■ 10"^ (see Fig. H]). This is not surprising because the total entropy should 
be proportional to the total multiplicity of secondary mesons which is approximately same 
for all good fits. As one can see from the last column of Table UTTt the corresponding E/S- 
ratios fall into a narrow interval 0.86 — 0.89 GeV. This observation is similar to the result 
of Ref. 451] that the observed ratio of the rest-frame energy to the multiplicity of produced 
hadrons is constant as a function of the bombarding energy. 

To check the sensitivity to the parameters of the phase transition, we have also calcu- 
lated the pion and kaon rapidity distributions for the EOS-II. It is found that with the 
same initial energy profiles as for the EOS-I, it is not possible to reproduce the observed 
spectra at any freeze-out temperature. In particular, the predicted kaon yield is strongly 



overestimated [6l|] at 100 MeV < Tp < Th = 190 MeV. Nevertheless, the BRAHMS data can 
be well reproduced with the EOS-II too, when taking smaller initial energy densities as com- 
pared with the EOS-I. Fits of approximately same quality are obtained for eo — 5 GeV/fm^ . 
As before, in choosing the initial conditions we apply the constraint fl23|) for the total energy 
of produced particles. Similarly to the case of the EOS-I, the data are better reproduced 
for initial profiles with small ?7o ^ 1 • 

Figures [HHTS] show the model results for pion and kaon rapidity distributions obtained 
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for the EOS-I and EOS-II. These results correspond to Gaussian initial profiles with 770 = C . 



For both EOSs we choose the parameter eo to obtain the best fit of the BRAHMS data|62|. 
Although the overall fits are very similar for both EOSs, the rapidity spectra obtained 
with the EOS-II are slightly broader than those with the EOS-I. In the same figures we 
demonstrate sensitivity to the choice of the freeze-out temperature. The best fits of the 
pion spectrum are achieved with Tp ~ 130 MeV (see Fig. [T^ . On the other hand, the kaon 
spectrum can be well reproduced only by assuming that kaons decouple at the very beginning 
of the hadronic stage, i.e. at Tp ~ T/^ = 165 (190) MeV for EOS-I (II). The contribution of 
resonance decays turns out to be rather significant, especially in the central rapidity region, 
where it amounts to about 35% (45%) of the total pion (kaon) yield. 
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FIG. 16: Rapidity distribution of vr '^— mesons in central Au+Au collisions at -y/sjvjv = 

200 GeV 

calculated for initial conditions with parameters eo = 50GeV/fm^, rjQ = {), a = 0.96. The solid 
line corresponds to the EOS-III. The dashed line is calculated for the hadronic EOS P = e/3. 



Experimental data are taken from Ref. 



According to Fig. [TU larger yields of secondary pions are predicted for smaller freeze- 
out temperatures. A much weaker sensitivity to Tp is found for kaons (see Fig. [T5l) . This 
difference can be explained by the large difference between the pion and kaon masses. Indeed, 
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in the case of direct pions, a good approximation at Tp > 100 MeV is to replace the transverse 
mass rriT in Eq. fl29l) by the pion transverse momentum pt- Neglecting the term with Tp in 
the r.h.s. of this equation, one can show that the rapidity distribution of pions at ?/ = is 
proportional to ^ = Tp{r]) cosh ?7-T|_/ cosh^ ^f(^) integrated over all i] . For a rough estimate, 
one can use the Bjorken relations j2| Yp = rj, spTp = sqTo , where sp is the entropy density 
at T = Tp . Using Eq. (ITT!) one gets Tp oc Sp^ oc Tp^^^" and therefore, ^ oc t|. , This 



shows that for < 1/3 the pion yield grows with decreasing Tp. Qualitatively, one can 
say that at low enough ch the increase of the spatial volume at freeze-out compensates for 
the decrease of the pion occupation numbers at smaller Tp. This effect is somewhat reduced 
because of a decreasing resonance contribution at smaller temperatures. It is obvious that 
for kaons this effect should be much weaker due to the presence of the activation exponent 
exp {—niK /Tp) . In fact, a numerical calculation for the same EOS and initial state shows 
that the kaon yield changes nonmonotonically: first it slightly increases when temperature 
goes down but then it starts to decrease at lower Tp . 
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FIG. 17: Same as Fig. (TUJ but for i('+-meson rapidity distributions. 



To investigate sensitivity of particle spectra to the presence of the phase transition, we 
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have performed calculations with purely hadronic EOSs. In this case we use the same initial 
conditions as before, but apply Eq. f|TT]) for all stages of the reaction, including the high 
density phase. Our analysis shows that for soft hadronic EOS with ajj < 0.2 it is possible to 
reproduce the observed pion and kaon data with approximately the same fit quality as in the 
calculations with the quark-gluon phase transition. Furthermore, the corresponding freeze- 
out temperatures do not change significantly. However, we could not achieve satisfactory fits 
for the "hard" hadronic EOS with cj^ ^ 1/3. Figures [T614T71 show particle spectra calculated 
for the same EOSs and initial conditions as in Fig. [TTl One can see that the hadronic EOS 
with cf^ = 1/3 leads to too wide rapidity distributions for both pions and kaons. The reason 
is that the higher pressure in the hadronic EOS gives a stronger push to the matter in 
forward and backward directions. From these findings we conclude that a certain degree of 
softening of the EOS is required to reproduce the pion and kaon rapidity distributions. Our 
analysis shows that reasonable fits of data for EOSs with ch ^ 0.2 (with and without the 
phase transition) can be achieved only when freeze-out temperatures for kaons are smaller 
than for pions. We find this situation unphysical, because kaons have smaller rescattering 
cross sections and, therefore, there is no reason for delayed kaon emission. 
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FIG. 18: Same as Fig. [T^ but for antiproton rapidity distributions. All results are obtained 
assuming fip = 0. Experimental data are taken from Ref. 28l |. 



It turned out that using the parameter sets from Table one can also reproduce rea- 
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sonably well the antiproton rapidity spectra measured by the BRAHMS Collaboration |28 |. 
Figure [18] shows the antiproton rapidity distributions, calculated for the parameter set A. 
In this case we explicitly take into account the contribution of the A(1232) irp decays, 
ignoring the width of A-isobar. Contributions of higher antibaryon resonances are taken 
into account in a similar way as for pions and kaons. The resonance contribution reaches 
about 55% at Tp = 165 MeV. One should consider these results as an upper bound for 
the antiproton yield. A more realistic model should include the effect of nonzero baryon 
chemical potentials which will certainly reduce the antibaryon yield. The thermal model 



analysis of RHIC data, performed in Refs. 



1471 ] ■ gives rather low values for the baryon 



chemical potentials, /if ~ 30 MeV, at midrapidity. This will suppress the antiproton yield 
by a factor ~ exp {—hf/Tf) ~ 0.8 . 

C. Rapidity distribution of total energy 




FIG. 19: Rapidity distribution of the total energy of secondary hadrons in central Au+Au collisions 
at ^Sjvjv = 200 GeV. Shown are results for the EOS-I and parameter set A. 
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In order to check the energy balance in the considered reaction, we have calculated 
additionally the rapidity distribution of the total energy of secondary particles, dE/dy. In 
the zero-width approximation this quantity can be written as 

f = cosh.X: Jd'PT^£^^mj+p^, (34) 

where the sum goes over all hadronic species. In calculating dE/dy we take into account 
not only direct pions (i = it) and kaons {i = K,K), but also 14 lowest mass mesons with 



widths r < 150 MeV 



M = r],p, u,ri', fo, ^o, 0, h, A, a2, v (1295), K*, K,, K; . (35) 

We also include 12 strange and nonstrange baryon-antibaryon (BB) species: 

B = N, Ar(1520), Ar(1535), A, A(1520), A(1620), A, A(1405), A(1520), S, S(1385), S . (36) 

In calculating contributions of different species to dE/dy in the r.h.s. of Eq. ( IMI) we 
use the freeze-out temperatures obtained by fitting the pion (Tp), kaon (Tp) and antipro- 
ton (Tp) rapidity distributions. We choose these temperatures to find, respectively, the con- 
tributions of free pions [i = n), heavy mesons [i = M , where M it) and baryon-antibaryon 
pairs (i = B,B). In particular, the results shown in Fig. [19] correspond to Tp = 130 MeV, 
Tp = Tp = 165 MeV. By integrating dE/dy, we have determined Ei and E^, the total 
energies of secondaries within the rapidity intervals |?/| < 1 and \y\ < 3, respectively. The 
BRAHMS Collaboration estimated Ei^^ from the rapidity distributions of charged pions, 
kaons, protons and antiprotons in most central Au+Au collisions at ^/s^ = 200 GeV. The 



values El ~ 1.5 TeV, E^ ~ 9TeV have been reported in Ref. 48] . From Table UTTl one can 
see that these values are well reproduced by the model. 

Based on the above analysis we conclude that within the hydrodynamical model the 
BRAHMS data can be well described with the EOS-I and EOS-II and the parameters of the 
initial state [tq = 1 fm/c) cr ^ 1.5 — 2, "^^o ^ 1 • The maximal initial energy density, eo, is 
sensitive to the critical temperature of the phase transition. For the EOS-I (T^ ~ 167 MeV) 
we get the estimate eo ^ 9 ± 1 GeV/fm^ while for the EOS-II (Tc ~ 192 MeV) the required 
values of eg are lower by about a factor of two. 
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D. Comparison of Landau like and Bj or ken like initial conditions 

As discussed above, the initial profiles which give the best fits of particle rapidity distribu- 
tions are intermediate between the Landau and Bjorken limits. Nevertheless it is interesting 
to compare the model predictions for the Landau-like and Bjorken-like initial states. In 
fact, the Landau initial conditions correspond to points situated near the origin of the (J — t]q 
plane in Fig. HI On the other hand, the Bjorken-like states correspond to table-like states 
with small a and large tjq i.e. to points near the vertical axis in Fig. H], away from the origin. 

Let us first consider the limiting case of the Landau model. In its simplest version, for 



T=1 fm/c 




FIG. 20: Temperature profiles at different proper times r calculated for the EOS-111. Left (right) 
panel corresponds to the Gaussian (table-like) initial conditions with parameters given in the 
figures. 



the EOSs with Cg = const, this model predicts a Gaussian rapidity distribution of secondary 
pions. Originally, the Landau model was formulated for proton-proton collisions. Below we 
assume that pion rapidity distributions in pp-reactions and central collisions of equal nuclei 
have the same shapes at the same cm. energy Then the following approximate 

formulae can be derived 



dy 



s 

^ '^iVJV 



0"La 



exp 



(37) 
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where the normahzation constant C-,^ does not depend on the bombarding energy and 



a 



r2 

Lan 



In 



(38) 



2(1 + c2) ' 31 -cf" 2mAr 

For = 1/3 and ^/s^ = 200 GeV one obtains a = 1/4 and aLan — 2.16 , the values often 
quoted in the hterature 0, [s^. On the other hand, for = 0.15 , which is preferable from 



the viewpoint of our previous analysis (see Sect. IIIIBI) . the predicted width is only 1.38 
i.e. significantly smaller than observed by the BRAHMS Collaboration. Note, that the 
analytical formulae fl371) - fl38l) are obtained for the case of constant sound velocity and can 
not be applied for EOSs with a first-order pase transition. In this connection we think that 



attempts 



50| to extract a unique sound velocity from pion spectra observed at a specific 



bombarding energy are unjustified. 
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FIG. 21: Same as Fig. [20l but for collective rapidity profiles. 



Within the Landau model one assumes that a thin disk-like fireball is created at some 
time t = to after beginning of the nuclear collision | 



It is postulated that this fireball is 
initially at rest {Y = 0) in the cm. frame. It has the transverse cross section S± = nR^ 
and the longitudinal size A2; ~ 2/2/7 • Here R is the geometrical radius of colliding nuclei, 
7 ~ \/s^/2mN is the cm. Lorentz-f actor. The total energy, E , is normally assumed to be 
E = ^{■y — 1) A^part^^Af where Apart is the number of participating nucleons and n ~ 0.5 — 0.7 
is the so-called inelasticity coefficient. At the cm. energy ^/s^ = 200 GeV the Lorentz- 
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factor 7 ~ 100 and, therefore, A2; < 0.1 fm. If to is not too small, the initial fireball 
occupies a region of small space-time rapidities, |?7| < Ar/ ^ Az/to, and one can approxi- 
mately apply the parametrization of initial state ([H]) with 770 = and cr ~ 2i?/(7to)- Note, 
that at small |?7| the difference between the initial flow proflles Y = rj and F = is negli- 
gible. Using further Eq. fl24p . one can estimate the maximal energy density of the fireball 
as eo ~ (v2vr^-R^o"to)~^ • On the basis of these arguments we conclude that the Landau 
model can be approximately simulated within our approach by taking the Gaussian initial 
profile with small a oc 7"^ and large eo oc 7^ . 
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FIG. 22: Isotherms in the rj — t plane calculated for the initial parameters eo = 782 GeV/fm^, 
r/o = , 0" = 0.1. Left panel corresponds to the EOS-III. Right panel shows results for the hadronic 
EOS P = e/3. Shaded region indicates the mixed phase. Dots near the origin correspond to the 
maximal initial temperature. 

For illustration, in this section we present the numerical results for the parameters of 
the initial state rjQ = 0, a = 0.1 and to — tq = Ifm/c. Again, the total energy is fixed at 
the value (12^ . This gives the estimate eo — 782GeV/fm'^. Figures [2DH2I] show profiles of 
the temperature and fiow rapidity at different proper times r ^ tq . For comparison with 
the Bjorken-like scenario, we also present results for a "broad", table-like initial state with 
cr = 0, ?7o = 2.3, eo = 20GeV/fm^ (this state has the same total energy and entropy). The 
difference between the fluid dynamics predicted for the Landau-like and Bjorken-like initial 
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states is clearly visible at r < lOfm/c and I?]] > 2 . It is interesting that asymptotically, for 
r > lOOfm/c, the "Landau" profiles at \ri\ < 2 behave similar to Bjorken's scaling solution, 
y ~ r], T ~ const . On the other hand, the two solutions differ significantly at large r] , 
where the Bjorken-like initial state leas to a very strong acceleration of the fluid near its 
outer edge. 

Figure [221 shows contours of constant temperature for the Landau-like initial state. These 
contours are calculated for two EOSs, with and without the phase transition. Contrary to 
the case of a broader initial energy density profile presented in Fig. [121 iiow the contours have 
a "pocket-like" structure with back-bending portions of isotherms. One can see a moderate 
sensitivity of isotherms to the presence of a phase transition. As discussed in Sect. IIII A[ 
this sensitivity becomes even weaker for softer EOSs. 

Now we present the rapidity distributions of pions and kaons calculated with the Landau- 
like initial conditions for several EOSs. Our analysis shows that it is not possible to repro- 
duce the pion spectra with EOS-I and EOS-II. The corresponding pion yields are too high 
at any Tp < Th ■ Both the pion and kaon yields become acceptable only for EOSs with 
cjj > 0.2. This is demonstrated in Figs. [23] - [24l which show the pion and kaon spectra for 
two EOSs with = 1/5 and 1/3. The first EOS corresponds to the EOS-I, but with 
cjj = 0.2 (this is in fact the EOS LII12 from Ref. [13l])- In Fig. [23[one can see that the pion 
spectrum calculated with this EOS is somewhat broader as compared with the experimental 
data. On the other hand. Fig. [24[ shows that for the EOS-III {cjj = 1/3) the predicted 
kaon spectrum is noticeably broader than observed experimentally. Therefore, we conclude 
that the BRAHMS data disfavor the Landau-like initial conditions at RHIC bombarding 
energies. 

Now let us consider in more details the Bjorken-like initial conditions corresponding to 
the table-like profiles {a = 0) of energy density 6J]. In this case the energy density at 
r = To is parametrized as eo6(?7o — jr^l) . At fixed ?7o the parameter eo can be determined 
from the total energy constraint ( l23l) . The explicit expression for rjo as a function of eo 
is given by Eq. ( I26i) with a = 0, f3 = 1 . Note, that strictly speaking, the Bjorken limit 
corresponds to infinite total energy and rjo ^ oo . Our calculations show that deviations 
from the Bjorken scaling increase with decreasing tjq . 

Figure [25] shows isotherms for a table-like initial state corresponding to the point D in 
Fig. [H Similarly to the case of Landau-like initial conditions, these isotherms have parts 
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FIG. 23: Rapidity distribution of 7r~'~— mesons in central Au+Au collisions at y^Sjvjv — 

200 GeV. 

Shown are results of hydro dynamical calculations for initial state with parameters 



eo = 782 GeV/fm'^, rjQ = 0, a = 0.1 . Experimental data are taken from Ref. 



29| 



with inverse slope, but they are much smaller and appear only at large rj . The calculations 
show that these parts contribute only to the tails of rapidity distributions. 

The pion and kaon rapidity distributions for Bjorken-like initial states are given in 
Figs. [271 Figure [21] shows the pion rapidity distributions calculated for several values 
of the parameter eo and for two EOSs with the phase transition. In all cases the freeze-out 
temperature is chosen by the best (when possible) fit of BRAHMS data. In accordance 
with the discussion in Sect. IIIIB[ the pion yield, calculated for fixed initial state, decreases 



(increases) with raising Tp for the EOS-I (III). 

The results for kaon spectra are given in Fig. [271 One can see that at eg ^ 10 GeV/fm^ the 
produced kaon yield is overestimated for the EOS-I at any freeze-out temperatures below 
Th = 165 MeV (we do not consider unrealistically small temperatures, Tp < 50MeV). On 
the other hand, at lower eo both pion and kaon spectra become too fiat at midrapidity. 
As seen from the right panels of Figs. [2614271 in the case of the EOS-III it is possible to 



34 



60 



40 - 



T3 



20 - 



T 



T 



T 



T 



T 



LH12, Tf=120 MeV 

EOS-III,Tf=147 MeV 



BRAHMS data: 




FIG. 24: Same as Fig. [23l but for ii'"'" -meson rapidity distributions. 
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FIG. 25: Same as Fig. [121 but for the initial profile with parameters eo = 7GeV/fm^ , a = 0, 
r?o = 3.34 . 



reproduce simultaneously pion and kaon spectra at cq ~ 20GeV/fm'^. However, for such 
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FIG. 26: Rapidity distribution of vr"'" -mesons in central Au+Au collisions at ^/sJ^ = 200 GeV. 
Solid, dashed and dashed-dotted curves are calculated for a = and different values of eo (in 
GeV/fm^), r/o and Tp (in MeV). The left (right) panel corresponds to the EOS-I (III). Experimental 
data are taken from Ref. [2^. 
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FIG. 27: Same as Fig. [26l but for ii'^'-meson rapidity distributions. 



eo the freeze-out temperatures for kaons are considerably lower than for pions. This means 
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that kaons are emitted later than pions, which is apparently unrealistic. 

On the basis of this analysis we conclude that both the Landau-like and the Bjorken-like 
initial conditions lead to unsatisfactory results for all considered EOSs. 



IV. SUMMARY AND DISCUSSION 

In this paper we have generalized Bjorken's scaling hydrodynamics for finite-size pro- 
files of energy density in pseudorapidity space. The hydrodynamical equations were solved 
numerically in r — r/ coordinates starting from the initial time tq = 1 fm/c until the freeze- 
out stage. The sensitivity of the final particle distributions to the initial conditions, the 
freeze-out temperature and the EOS has been investigated. A comparison of vr, K, p rapid- 
ity spectra with the BRAHMS data for central Au+Au collisions at ^/s^ = 200 GeV has 
been made. The best agreement with these data is obtained for initial states with nearly 
Gaussian profiles of the energy density. In choosing the initial conditions we impose the 



constraint on the total energy of produced particles known from experimental data 28|. It 
is found that the maximum energy density of the initial state, eo, is sensitive to the param- 
eters of a possible deconfinement phase transition. The BRAHMS data are well reproduced 
with eo of about 10 (5) GeV/fm^ for the critical temperature Tc ~ 165 (190) MeV. The only 
unsatisfactory aspect of these calculations is the prediction of very long freeze-out times, 
~ 50 fm/c, for pions. 

We would like to comment on several points. First, all calculations in this paper are made 
for tq = 1 fm/c. Of course, one can start the hydrodynamical evolution from an earlier time, 
i.e. assuming smaller tq . In this case one should choose accordingly higher initial energy 
densities. But tq cannot be taken too small, since at very early times the energy is most 



likely stored in strong chromofields 



51| . The quark-gluon plasma is produced as a result of 



the decay of these fields. Estimates show that the characteristic decay times are in the range 



0.3 — 1.0 fm/c. At earlier times the system will contain both the fields as we^ 
partons, and the evolution equations will be more complicated, see e.g. Ref. 



1 as produced 



52|. 



It is obvious that the Cooper-Frye scenario of the freeze-out process, applied in this 



paper is too simplified. This was demonstrated e.g. in Ref. 531]. One should also bear 
in mind that the freeze-out temperatures obtained in our model will be modified by the 
effects of transverse expansion and chemical nonequilibrium. Attempts to construct a more 
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realistic description of the freeze-out have been performed in Refs. |ld . Il3| . where the 
transport model was applied to describe evolution of the hadronic phase. In this approach 
the solution of fluid-dynamical equations is used to obtain initial conditions for hadronic 
cascade calculations at later stages of matter expansion. We are planning to study these 
problems in the future. 



APPENDIX A: THE FLUID DYNAMICAL EQUATIONS IN THE CONSERVA- 
TIVE FORM 



In order to apply the standard flux-corrected transport algorithm 3J] to solve numerically 
the equations of fluid dynamics it is necessary to represent them in the conservative form. 
The following set of equations can be written within the 1+1 dimensional model in the 
cartesian coordinates t — z (for completeness we consider the general case of nonzero baryon 
charge) 

dtu + d,{uv, + f) + g = 0, (Al) 

where z and Vz = U^/U^ are the spacial coordinate and the fluid 3-velocity. The multi- 
component vector u includes the densities of baryon charge {ui = nU^), energy {u2 = T°°) 
and momentum (^3 = T°^). In the usual space-time representation g vanishes and /i = 0, 

f2 = PVz, h = P. 

To rewrite these equations in the light-cone variables, one can perform the Lorentz- 
transformation of the baryon current = nU^ and the energy-momentum tensor T^"^ to 
the frame moving with the rapidity r] = tanh~^(z/t) (this reference frame is marked below 
by tilde). In particular, J° = ncoshY, = nsinhY and 

f °° = e cosh^ Y + P sinh^ Y , (A2) 
f = (e + P) cosh Y sinh Y , (A3) 
f'^ = e sinh^ Y + P cosh^ Y , (A4) 

where Y = Y — r] . 

From Eqs. ([I])-([2]) of the main text we have 



drJ^ + d^J' + J^T = 0, (A5) 
<9,f°° + <9^f°^ + (f°° + f"')/r = 0, (A6) 
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drT'^' + d^T'' + 2T°7r = , (A7) 

where dx = rdr] . These equations are equivalent to Eqs. ([5])-(j7]) of Sect. Ill A[ One 
can see that Eqs. flA5p - flA7p have the form flAip with the replacement t —>■ r, z ^ tt], 
Vz — > tanh(Y — rf). 

In the baryon-free case J° = = and Eq. ( 1A5I) is trivially satisfied. The "phoeni- 
cal" version of SHASTA algorithm 3J] is used to solve numerically Eqs. (1A6l) - (1A7ll in the 
rectangle \rj\ < r^max, To < t < Tmax of the t] — r plane. Typically we choose r^max = 10, 
Tmax = 200fm/c. The grid size Arj = 1.25 ■ 10~^ is used in most calculations. For each time 
step the increment At is found from the requirement At = min(0.1rAr7, O.lfm/c) . The 
latter is motivated by the Courant condition for the stability of fluid-dynamical equations 
in a finite difference approximation. Within the light cone representation this conditions 
takes the form At < Ax/cg = tAti/cs , where < c is the sound velocity. 



APPENDIX B: PARTICLE SPECTRA FROM RESONANCE DECAY 

One can simplify the expressions fl30l) for spectra of particles produced in resonance 
decays, having in mind that degeneracy effects in resonance phase-space densities should 
be small at small and realistic temperatures Tp > 100 MeV. Therefore, when calculating 
invariant momentum distributions of resonances, one can approximate the Bose or Fermi 
functions in Eq. (127|) by the corresponding Boltzmann-Jiittner expression: 

where qr is the degeneracy factor of the resonance R . Presence of delta-function allows us 
to easily remove integration over the angle between and p in Eq. ( l30l) . Using further 
the Lorentz-invariance of ( 130|) . it is possible to make analytic integration over the resonance 
energy in the fluid rest frame. The result can be written as 

oo 

d^NR^i d^Nn^i QRbR^^ f tur \ f ^ [a u , orrMN t-no\ 

The Lorentz-scalars A and B in the r.h.s. are determined as follows 



A = ^ [E'iE. + Tp) - mREo.]e(^'P ' ^-)/^^ - -> , (B3) 



Tf 



5 = -| [E'ttirEo. - m1{E^ + Tf)] e(/^^ " ^-)/^^ - {E^ ^ E+) , (B4) 
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where E' = ^/mf + p''^ = pUp = mTCOsh.{y — Yp) and E± = mji{EQiE' ± qop')/^^ ■ 

In most cases, when calculating contributions of resonance decays we apply the zero- 
width approximation, i.e. substitute w{m^) = 6{mji — mji) , where m/j is the central mass 
value [gS] presented in Ref. 42 1 for the resonance R. However, we take into account the 
nonzero width of p-mesons. Following Ref. 40|, we apply the parametrization 



wimR) 



(B5) 



Here the constant C is taken from the normalization condition, J dmiiw{mii) = 1 , and 



(B6) 



where = ^/vp^^IAl — . The numerical values 40] 



mR = 770 MeV, A = 159 MeV, ^ = 0.5 , 



(B7) 



are used to calculate feeding from p-meson decays. 
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